Absence of Asymptotic Freedom in Doped Mott Insulators: Breakdown of Strong 

Coupling Expansions 
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We show that doped Mott insulators, such as the copper-oxide superconductors, are asymptoti- 
cally slaved in that the quasiparticle weight, Z, near half-filling depends critically on the existence 
of the high energy scale set by the upper Hubbard band. In particular, near half filling, the following 
dichotomy arises: Z / when the high energy scale is integrated out but Z — in the thermody- 
namic limit when it is retained. Slavery to the high energy scale arises from quantum interference 
between electronic excitations across the Mott gap. Broad spectral features seen in photoemission 
in the normal state of the cuprates are argued to arise from high energy slavery. 



Vast progress in theoretical solid state physics has been 
made by constructing models which describe the low- 
energy properties of solids. Essential to the success of 
this program is the separability of the high and low en- 
ergy degrees of freedom. Should this separability hold, 
then the high energy scales can be integrated out yield- 
ing an effective Hamiltonian which describes the relevant 
low-energy or long-wavelength physics. Notable successes 
include the prediction that dilute magnetic moments are 
quenched at low temperatures in non-magnetic metals 
(the Kondo problem) [l| and the Landau Fermi liquid 
theory|2( of the normal state of metals. 

In the context of high-temperature superconductivity 
in the copper oxide materials, a similar approach has 
been adoptedH H IE S H ES E3 The high energy 
scale in the cuprates corresponds to the energy cost to 
doubly occupy the same 'lattice site' in the copper-oxide 
plane. At half-filling, an insulating state (the Mott in- 
sulator) obtains when the double occupancy energy cost 
(U) vastly exceeds the nearest-neighbour hopping energy, 
t. In the cuprates, U ~ lOt. This corresponds to the 
strong coupling regime. In this regime, a requirement for 
any low-energy theory is that all hopping processes pre- 
serve the number of doubly occupied sites The resul- 
tant Hamiltonian can then be projected onto the singly 
occupied subspace or lower Hubbard band and then stud- 
ied accordingly [J 00 00, EH E3- Through second or- 
der in t/U, the spin-spin Hamiltonian (or t — J model) 
that obtains has been used widely as an effective model 
for the cuprates 0,0, 0,0,0- In this Letter, we show that 
this procedure fails for the high-temperature supercon- 
ductors because in the strong-coupling regime, the value 
of the quasiparticle weight, in the vicinity of half-filling, 
depends crucially on the presence of the upper Hubbard 
band. Namely models in which it is integrated out yield 
well-defined quasiparticles, whereas those in which it is 



not yield a vanishing weight in the thermodynamic limit. 
At work here is a non-perturbative quantum interference 
effect between excitations that live at energy scales that 
span the Mott gap. 

Our conclusion that the physics of projected models 
can be different from that of the Hubbard model suggests 
that the high and low energy degrees of freedom are cou- 
pled. Indeed, numerous experiments support this view. 
For example, optical conductivity measurements [l^. E5j | 
as well as oxygen K-edge photoemissionE3 indicate that 
there is a massive re-shuffling of spectral weight from an 
energy scale as high as 2eV above the Fermi energy in 
both hole and electron-doped cuprates as a function of 
doping such that the low-energy spectral weight (LESW) 
increases at the expense of the high-energy (> 2eV) spec- 
tral weight. Similar results are also seen in angle resolved 
photoemission The non-trivial sum rule 0,0] that 
has emerged from oxygen K-edge x-ray studies is that at 
least one single particle state (per doped hole) is lost at 
high energies and transferred to low energies such that 
the LESW increases faster than 2x, where x is the hole- 
doping level. The t — J model does not preserve E3 
this sum rule as the LESW is exactly 2x in this trun- 
cated scheme. LESW in excess of 2x is purely dynamical 
and arises from virtual excitations to the upper Hub- 
bard band. Such spectral weight transfer indicates that 
some low-energy degrees of freedom in the normal state of 
the cuprates are derived from a high energy scale. The 
strong electron correlations that give rise to a mixing 
between the low and high energy degrees of freedom in 
doped Mott insulators we termed Mottness\V^ . A further 
surprise 0, 0, |2(| is that Mottness persists even when 
su perc onductivity obtains. For example, Rubhaussen, et. 
al.[18j have shown that changes in the optical conductiv- 
ity occur at energies 3eV (roughly 100A, A the maximum 
superconducting gap) away from the Fermi energy at T c , 
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and Bolegraf, et. al.Jjj have seen an acceleration in the 
depletion of the high energy spectral weight accompa- 
nied with a compensating increase in the low-energy spec- 
tral weight at and below the sup erconducting transition. 
Similarly, Bontcmps, et. al. 2(j observed that in under- 
doped (but not overdoped) cuprates, the Glover-Ferrel- 
Tinkham sum rule is violated and the optical conduc- 
tivity must be integrated to ~ 100A, as opposed to 4 A 
in conventional superconductors, for the spectral weight 
lost upon condensation into the superconducting state to 
be recovered. 

Nonetheless, the experiments dllE E3, EE E 13 E3 
which demonstrate that all energy scales are mixed in the 
cuprates have had virtually no theoretical impact. What 
seems to be missing is an explicit demonstration that a 
key physical quantity differs once the high energy scale is 
integrated out. We propose that a crucial quantity that 
captures the difference between projected models and the 
Hubbard model is the quasiparticle weight. 

To proceed, we first show that all problems regarding 
the formulation of projected models to high order can 
be overcome. Hence, should a failure arise, it does not 
reside in the formulation. Consider the simplest model 
for a doped Mott insulator. In the Hubbard model, 

H = T + V = -tJ2 el*. + E/ £ ni t mi (1) 

(ij) 1 

electrons acquire kinetic energy by hopping among neigh- 
bouring sites, (ij) and experience an on-site repulsive in- 
teraction, U. It is expedient to break the kinetic energy 
into three terms, 

To = -tE(iv+4^). 

(ij) 

T l = -t^ViaZjr, 

(ij) 

T-i = (2) 

(ij) 

which are eigenoperators of the interaction and as a re- 
sult, obey the commutator, [V, T m ] — mUT m . The op- 
erator T m increases the double occupancy by m. The 
operators r] irT = c^n^^ and £, CT = c icr (l - 7ii_ CT ) anni- 
hilate electrons on doubly and singly-occupied sites, re- 
spectively. Note, Ci a = rjiu + £ i(T . Successful removal of 
double occupancy implies that r\i„ and can be decou- 
pled. For simplicity, we will set U = 1. 

To show how such excitations enter the projected 
Hamiltonian schemes, wc review the two standard per- 
turbative approaches used in this context. In the first ap- 
proach to removing double occupancy, we use a similarity 
transformation, S, such that the transformed Hamilto- 
nian 

H = e s He- s . (3) 



does not contain hops between sites with differing num- 
bers of doubly occupied sites. As this procedure is well- 
described[ij|2l|,|22|,|23| in the literature, we will be brief. 
Our Hamiltonian initially is H = V + T + T_i + Z\. 
The last two terms in the Hamiltonian do not con- 
serve the number of doubly occupied sites and hence 
must be eliminated. In the standard implementation, 
the similarity transformation]^ is chosen such that 
T-i + Ti + [V, S* (1) ] = 0. The transformation that ac- 
complishes this is S 1 -- 1 ^ = Ti — T-i. At each order, the 
similarity transformation must be modified accordingly. 
To obtain the effective Hamiltonian in the singly occu- 
pied subspace, we perform the projection P HP which 
removes all terms in which T_i appears first as a conse- 
quence of T_i\Q) = 0, where |0) is any state in the lower 
Hubbard band. Through fourth order, we write the ef- 
fective Hamiltonian j2^|, H e g — iTcd + H 1Ile< i, as a sum 
of irreducible, 

ffirrcd = 70-^(1 + T + T 2 )T 1 +T! 1 T 1 2 (4) 

and reducible, 

ffred = {T - Tl ^y^} + {To, T_iT)Ti} + (T_!T!) 2 

terms where {a, b} = ab + ba. Third and higher order 
terms are non-zero away from half-filling in bi-partite 
lattices. Since each term preserves the number of doubly 
occupied sites, the sum of the indices on each product 
of T„'s vanishes. In the irreducible terms, all the inter- 
mediate states contain at least one doubly occupied site. 
The energy denominators ( the 1/U factors) arise from 
this energy difference. The reducible terms are prod- 
ucts of irreducible ones and hence they contain hopping 
processes that do not originate from excitation to the 
doubly occupied subspace, T_iTiTo nor terminate once 
an electron is returned to the singly occupied subspace, 
for example, (T^T) 2 . All such processes can be viewed 
as arising from a rotation|24| of the eigenstates in the 
low-energy sector. Such a rotation arises naturally in 
this context, since a unitary transformation preserves or- 
thogonality. While it might be anticipated that the true 
effective low-energy Hamiltonian should be independent 
of such a rotation and hence irreducible with respect to 
the target manifold, such is not the case here. In fact, 
the super-extensive parts of iT rr cd cancel those of -ff re dj 
giving rise to a linked expansion for the energy. The first 
two terms in -ffi rre d yield the t — J model in addition to 
the three-site hopping which describes the motion of a 
hole in a spin background. 

The relationship between the change of basis and size 
consistency is further illustrated using Brillouin-Wigncr 
(BW) perturbation theory. Let P be the projector for 
the lowest degeneracy subspace and Q = 1 — P, the or- 
thogonal complement. Because V and To do not change 
the number of doubly occupied sites, [Q, V + Tq] = 0. 
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Consider the Schrodinger equation, (E — V — T)\ip) = 0, 
where \ip) is the exact many-body eigenstate in the Hub- 
bard model and P\ip) — \ipo) yields the exact eigenstate 
in the zero-double occupancy sector. Multiplying the 
Schrodinger equation on the left by Q results in the for- 
mal expansion (E-V-T )Q\ip) = Q(T_i+Ti)|V> for the 
components of the eigenstates orthogonal to those in the 
zero double occupancy subspace.To compute the energy 
eigenvalue, we multiply the Schrodinger equation on the 
left by P to obtain E\i) ) = T \tpo) + PT^Q\ip). Succes- 
sively iterating this equation twice by using the equation 
for Q\^p) yields the self-consistent expansion 



E\ipa) = \T + PT- X 



E-V-To 

1 ii- 1 



E-V-Tn E-V-T, 



E-V -T 
Ti ) \iPo) (5) 



for the energy through fourth order. To obtain a more 
useful form for the energy eigenvalue, we expand the 
energy denominators and multiply on the left by (ipo\. 
What results is a second-order polynomial in E. The 
unique root that vanishes as t — > oo is given by 

E = (H imd ) - (To - T_iTi)(T_iTi) + {T f {T^T{) (6) 

through 0{t A ). The expectation value in Eq. (0 is per- 
formed with the exact eigenstates in the subspace with 
zero double occupancy. Typically in degenerate pertur- 
bation theory, a basis which lifts the degeneracy to first 
order is sufficient to evaluate all the higher order terms. 
In this case, this would correspond to using a basis that 
diagonalizes To. In such a basis, T-\T\ is not diag- 
onal and hence the unlinked part of (T_iT_iTiTi) = 
((T_iTi) 2 ) cx N 2 is not canceled. This problem is en- 
demic to the Hubbard model, because in traditional per- 
turbation theories there is no analogue of To which in- 
duces transition only in the target space. The correct 
scaling with N is accomplished by expanding \ipo) in pow- 
ers of t and collecting all unlinked terms order by order in 
t. If this is done, the unlinked part of the terms contain- 
ing a single factor To and (T_iTi) 2 in iJi rre d cancel the 
second term in Eq. (0 and the unlinked part of the Tq 
terms in i?i rr ed cancel the last term in Eq. © through 
fourth order. The unlinked parts of these terms of course 
have 0(t 5 ) contributions and higher. All of these terms 
can be shown to cancel by the order by order expansion 
of the eigenstates in the lowest energy sector. Conse- 
quently, perturbation theory and the canonical transfor- 
mation will be equivalent up to an arbitrary rotation in 
the target space. While such rotations affect -ff rc d not 
-ffirred, -ffred is crucial to the correct size dependence of 
the true effective Hamiltonian. 

While the subtleties in constructing projected effec- 
tive Hamiltonians can be overcome up to an arbitrary 
rotation in the target space, all such expansions rely on 



the partitioning of the electron into ^ and r\i a exci- 
tations. Consequently, the full electron spectral func- 
tion, A(k,co) = -lmFT{e(t - t')({c l(J (t),c] a {t')}))/n = 
A m + Ape + 2A V £, contains two diagonal terms corre- 
sponding to the upper and lower Hubbard bands, A m 
and A^, respectively and a cross term A V £ which rep- 
resents the degree to which the high and low energy de- 
grees of freedom are coupled. Here, FT represents the 
frequency and momentum Fourier transform. Shown in 
Fig. is an explicit calculation of the A v ^ term (in- 
tegrated over k) using the dynamical two-site method 
detailed previously 01- As m other cellular methods[25j|. 
the self-energy for the lattice is constructed from the re- 
solvents for the electronic states on a finite cluster using 
a self-consistent closure. In this case, a two-site cluster is 
used. As is evident from Fig. A,^ is distinctly non- 
zero and mirrors the overall single particle density of 
states with peaks at the upper and lower Hubbard bands. 
Three features are most relevant. First, at U — 8i, and 
at half-filling, the overall density of states (see Fig. 4 
of Phys. Rev. B, 64 235117/1-8 (2001)) has a maxi- 
mum value of 0.16, whereas the total weight arising from 
the cross term, 2A V ^, is .04 or 25%. Hence, this contri- 
bution cannot be ignored. Second, the cross term has 
both negative and positive contributions. This structure 
arises necessarily because the integral of A V £ over all fre- 
quency yields the equal time correlator ({£i<7, Via}) = 0' 
whose vanishing maintains the Pauli principle. This im- 
plies that A v £ is either zero, which it is not, or it must 
have both positive and negative parts, representing con- 
structive and destructive interference, respectively, be- 
tween different regions in energy space. Third, when U 
is increased by 50% to U = 12t as in the cuprates, the 
cross term does not decrease appreciably. Nonetheless, 
at half-filling, the contribution of the crossterm below 
and above the chemical potential sum to zero indepen- 
dently, indicating that the upper Hubbard band can be 
integrated out safely without sacrificing the Pauli princi- 
ple. 

Such is not the case, however, at finite doping. The 
lower panel in Fig. indicates that a pseudogap de- 
velops at the chemical potential, indicating an orthog- 
onality catastrophe and hence a vanishing of the quasi- 
particle weight. Because a pseudogap subtracts spectral 
weight at low energies and transfers it to the upper Hub- 
bard band, the sum rule which ensures the Pauli prinic- 
ple is satisfied only when A V £ is integrated over all en- 
ergy scales not simply up to the chemical potential (or 
some intermediate energy cutorr) as would be case in 
projected models. Hence, although the cross term can 
be obtained perturbatively from projected schemes by 
canonically transforming (as described previously) the 
electron operators, perturbation theory fails as the inte- 
grated cross term represents inherently non-perturbative 
physics, namely the Pauli principle. Symptomatic of this 
failure is the difference in the value of the single-hole 
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FIG. 1: Cross correlation or quantum interference be- 
tween the upper and lower Hubbard bands at half-filling, 
n — 1 and at n = .95 at T = O.lt. The dip at the chemical 
potential in the lower panel represents the pseudogap. The 
inset shows that this dip leads to a vanishing density of states 
at zero temperature and hence an orthogonality catastrophe. 



quasiparticle weight, Z = \(4'G\Ck cr \' l Pk,<j)\ 2 between the 
Hubbard and projected models. Here tpo and V'k.a- are 
the exact ground states for the half-filled and one-hole 
systems, respectively. In projected schemes, such as the 
t — J model, Z oc J/t as has been demonstrated both 
analytically and numerically [2^|. Consequently, a single 
hole is delocalized at T = in the t-J model. However, in 
the Hubbard model, adding a single hole leads to a "non- 
renormalizable" phase shift of each state in the first Bril- 
louin zone and hence an orthogonality catastrophe[4l l27j 
in the thermodynamic limit, Z oc L _l5 0, H3|- Conse- 
quently, Z for one-hole in the Hubbard model at half- 
filling does not appear to have a well-defined expansion 
in t/U. It is this breakdown that we term asymptotic 
slavery. This failure, applies strictly in the thermody- 
namic limit and should persist as long as a pseudogap 
is present which necessarily leads to a vanishing of the 
quasiparticle weight. Broad spectral features [l5| (that 
is, Z = 0), spectral weight transfer bole localization 
in the underdoped regime, as well as the color change 
seen in optical experiments [19j upon a transition to the 
superconducting state are all signatures of the quantum 
interference that is the root of asymptotic slavery. 

Of course, asymptotic slavery in doped Mott insula- 
tors stands in stark contrast to the perturbative physics 



present at short distances, that is asymptotic freedom, 
in quark matter |28l Ejj. Short of an exact construction 
of the quasiparticles, any realistic model of the cuprates 
must be solved on the energy scale U because double oc- 
cupancy does not necessarily mean high energy. That 
is, a Wilsonian renormalization group analysis fails as 
long as asymptotic slavery is present, namely as long 
as the pseudogap presists. The high energy scale does 
not simply renormalize the low-energy degrees of free- 
dom. In doped Mott systems, the Pauli principle appears 
as a sum rule over high and low energies. Methods 
which emphasize local non-perturbative physics or per- 
haps non-commutative field theories |3dl| (which display 
UV-IR mixing) have the ingredients to capture how Mot- 
tness conspires to yield asymptotic slavery. 
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